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Abstract. This paper derives completeness properties of the bispectrum for compact groups 
and their homogeneous spaces. The bispectrum is the Fourier transform of the triple correlation, just 
as the magnitude-squared spectrum is the Fourier transform of the autocorrelation. The bispectrum 
has been applied in time series analysis to measure non-Gaussianity and non-linearity. It has also 
been applied to provide orientation and position independent character recognition, as well as to 
analyze statistical properties of the cosmic microwave background radiation; in both cases, the data 
may be defined on a sphere. On the real line, it is known that the bispectrum is not only invariant 
under translation of the underlying function, but in many cases of interest, it is also complete, in 
that the function may be recovered uniquely up to a translation from its bispectrum. This paper 
extends the completeness theory of the bispectrum to compact groups and their homogeneous spaces, 
including the sphere. The main result, which depends on Tannaka-Krein duality theory, shows that 
every function whose Fourier coefficient matrices are always nonsingular is completely determined by 
its bispectrum, up to a single group action. Furthermore, algorithms are described for reconstructing 
functions defined on SU{2) and 50(3) from their bispectra. 
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1. Introduction. The triple correlation of a complex-valued function on the real 
line is the integral of that function multiplied by two independently-shifted copies of 
itself: 

/oo 
f*{x)f{x + si)f{x + S2)dx. 
-oo 

It is easily seen that the triple correlation does not change if the function is translated. 
The Fourier transform of triple correlation is the bispectrum. If the Fourier transform 
of / is denoted F = J-{f}, then the bispectrum is 

(1.2) ^3,/(«,«) = ^{a3,/(si,S2)} = F{u)F(v)F*{u + v) 

The triple correlation extends the concept of autocorrelation (denoted here 02,/) 
which correlates a function with a single shifted copy of itself, thereby enhancing the 
function's latent periodicities: 

/oo 
r{x)f[x + s)dx. 
-00 

The Fourier transform of the autocorrelation is A2J = J'{a2.f} = \F\'^ , which obvi- 
ously lacks phase information and therefore provides a limited analysis of a function's 
structure. In contrast, eq. (jl.2p shows that the bispectrum contains both magnitude 
and phase information, while still being invariant to translation. These properties 
suggest applications in invariant matching for pattern recognition. More importantly 
for matching, the bispectrum in many cases of interest is not only invariant but also 
provides a complete description of the function: the function may be reconstructed 
from it, up to a single unknown translation [llj . 

*This work was supported by NTU Research Support Grant M58020000. 

t School of Computer Engineering, Nanyang Technological University, Singapore. 



1 



2 



R. KAKARALA 



The bispectrum was perhaps first investigated by statisticians examining the cu- 
mulant structure of non-Gaussian random processes It is well known that the 
third order cumulant, of which the triple correlation is a sample, has zero expected 
value for a Gaussian process. Hence, the bispectrum is a tool for measuring non- 
Gaussianity. The bispectrum was also independently studied by physicists as a tool 
for spectroscopy. H. Gamo [6] in 1963 described an apparatus for measuring the triple 
correlation of a laser beam, and also showed how phase information can be completely 
recovered from the real part of the bispectrum — up to sign reversal and linear offset. 
Gamo's is perhaps the first completeness result along the lines of what is explored in 
detail in this paper. However, Gamo's method implicitly requires the Fourier trans- 
form to never be zero at any frequency. This requirement was relaxed, and the class of 
functions which are known to be completely identified by their bispectra was consider- 
ably expanded, by the study of Yellott and Iverson for example, every integrable 
function with compact support is completely determined, up to a translation, by its 
bispectrum. 

The statistical and physical applications described above are for data defined 
on Euclidean domains M". The bispectrum also finds applications on non- Euclidean 
domains such as the sphere S*^. For example, in astrophysics, the Cosmic Background 
Radiation (CBR) may be modelled as a function defined on a sphere. X. Luo [TB] 
calculates the bispectrum of spherical CBR functions and examines the properties 
the bispectral coefficients have under cosmological scenarios such as inflation and 
late-time phase transitions. 

The bispectrum's invariance and completeness have motivated researchers in pat- 
tern recognition to apply it for template matching and shape recognition. For exam- 
ple, R. Kondor [14, Ch. 8] shows how position and orientation independent optical 
character recognition can be accomplished by projecting the character from the plane 
on to the sphere, and subsequently using the bispectrum on the sphere for invariant 
matching. Kondor's results are discussed further in Section [5l 

Despite the interest in applying the bispectrum for non-Euclidean domains, little 
has been published about important properties such as completeness. The contribu- 
tion of this paper is to derive the completeness theory of the bispectrum for (non- 
commutative) compact groups and their homogeneous spaces. In order to construct 
the bispectrum on groups, we require concepts from harmonic analysis using group 
representation theory. Those concepts are presented in the next section. Then, a 
matrix form of the bispectrum which proves convenient for analysis is demonstrated. 
The matrix formulation allows a relatively simple criterion for completeness: it is 
shown that functions defined on compact groups that have nonsingular Fourier trans- 
form coefficients are completely determined by their bispectra. This result depends 
on the well-known Tannaka-Krein duality theory of compact group representations. 
The completeness result is extended to homogeneous spaces using the Iwahori-Sugiura 
duality theorem j9]. Reconstruction algorithms for functions defined on SU{2) and 
S0{3) in particular are described, expanding on material in a previous paper [13j . 

2. Preliminaries. Let us review the basic concepts of representation theory for 
compact groups. For more details, the reader may consult Chevalley ( 4 , Chapter 
VI). Let G be a compact group. A n- dimensional unitary representation of G is a 
continuous homomorphism D of G into the group U{n, C) of n x n unitary matrices. 
The following operations are all defined on representations: complex-conjugation, 
direct sum, and tensor product. The complex conjugate of any representation is simply 
that representation D* whose matrices are complex-conjugates of those of D. Let Di, 
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D2 be respectively a n-dimensional and a m-dimensional representation; their direct- 
sum Di © D2 is the (n + m)-diniensional representation which maps each element g 
to the block-diagonal matrix Di{g) © D2{g) with Di{g) in the upper left corner and 
D2{g) in the lower right corner. Similarly, the tensor product Di D2 of Di and 
D2 is nm-dimensional representation whose matrices are n x n blocks, each block of 
size m X m, where for each g the (i, j)-th block is the matrix D2{g) multiplied by the 
(i,j)-th coefficient of Di{g). We reserve the symbol 1 for the trivial representation 
that maps all of G into the number 1. 

We now define equivalence and reducibility of representations. Two representa- 
tions Di and D2 are equivalent if there exists a unitary matrix C such that Di{g) = 
CD2{g)C'^ for all g G, with f denoting the matrix conjugate transpose. A represen- 
tation D is reducible if there exists a matrix C and two representations Di, D2 such 
that D{g) = C [Di{g) © (.g)] G^ for all g G G; otherwise D is irreducible. Every 
representation D is the direct sum of irreducible unitary representations. 

The set Q of all equivalence classes of irreducible unitary representations of G is 
called the dual object of G. Intuitively, Q is the "frequency" domain of G. In general, 
Q is not a group, unless G is abelian. In what follows, we denote elements of Q, 
which are equivalence classes, by Greek letters such as a, /3, etc. For each a ^ G, 
let dim(a) denote the common dimension of the representations in a. Let {Da}aeg 
denote any set of representations that contain exactly one member in each equivalence 
class in G- We call any such set a selection. A given selection has two properties, 
both consequences of the classical result known as the Peter- Weyl theorem: (1) there 
are at most countably many representations in any selection, and in particular there 
are finitely many if, and only if, G is finite, i.e., G is finite if and only if G is finite; 
(2) the set of all matrix coefficients dPj^{-), where a E G and 1 < p,q < dim(a), from 
any selection forms an orthogonal basis for the Hilbert space L2{G). 

2.1. Duality theory. Our main tool for proving completeness results is the 
duality theorem due to T. Tannaka and M. Krein. There are several formulations 
of that result, of which Chevalley's [H pp 188-203] is the most convenient for our 
purposes. The formulation is as follows. Let 6(G) be the representative algebra of 
G, which is the algebra of complex-valued functions on G that is generated by the 
set of matrix coefficients of any selection {Da}ai£g. In fact, 6(G) is independent 
of the selection that is used. It is well-known that every function / G 6(G) may 
be expressed in exactly one way as a finite linear combination of the set of matrix 
coefficients from any given selection (01 pg 189]). The structure of 6(G) may be 
understood by considering algebra homomorphisms, which are maps oj : 6(G) — > C 
that are both linear and multiplicative, i.e., uj{cifi + C2/2) = ciuj{fi) + C2w(/2) and 
^{hf'i) = ^{fi)^{f2) for all scalars ci, C2, and functions /i, /2 in 6(G). The set 
of all algebra homomorphisms is denoted Q.{G). Clearly, for every g € G, the map 
ojgif) = f{g) is an algebra homomorphism. Note that ujg{f*) = ujg{f)*, i.e., ujg 
preserves complex-conjugation. In fact, the converse is also true, an identification 
that is essential to the duality between groups and their representations; see [H pg 
211] for details and proofs. 

Theorem 2.1 (Tannaka-Krein) . To every algebra homomorphism to G fl(G) 
that preserves complex- conjugation, i.e., Lo{f*) = uj{f)* , there corresponds a unique 
element g E G such that uj{f) = f{g) for all f G 6(G). 

We recast this duality theorem in a slightly different form (see also pDj pp 303- 
306]). Let {Da}aeg be any selection of irreducible representations, and let {U{a)}aeg 
be a corresponding sequence of unitary matrices, such that for each a the matrix U (a) 
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has the same dimension as the representation D^. We determine the necessary and 
sufficient conditions under which the latter sequence arises from the former by an 
element of Q{G), i.e., when U{a) = Lo{Da) for some fixed homomorphism u) G 57(G). 
Consider the tensor product D^^ ® Dg of any two representations in our selection; 
that representation is, in general, reducible, and we write its decomposition into 
irreducibles (taken from our selection) as follows: 

(2.1) D„®Ds^ C^s [D^, © ■ • ■ © D„J Cl^. 

The indices ai, . . ., au appearing on the right are unique up to permutation ([U 
pg 175]). Suppose now that there exists lo such that U{a) — uj{Da) for all a. By 
applying lu to both sides of eq. ()2.ip . and using the fact that co is both linear and 
multiplicative, we obtain that (writing U{a) for Lu{Da)) 

(2.2) U{a) ® U{d) = C,s [C/(ai) © • • • © [/(a^)] C^^ 

Equation (|2.2p is not only necessary, but also sufficient, as the following result shows. 

Theorem 2.2. Let {i^ajaee ^'^'^ {U{oL)}aeQ above. If, whenever eq. \2.1\) 

is true, we have that eq. i2.2\) is also true for same a, 6, and matrix CaS, then there 
exists a fixed g £ G such that U{a) = Da{g) for all a. 

Proof. To any sequence {U{a)}ag, there exists a unique linear map to : 0(G) C 
such that U{a) = uj{Da). To see this, note that the set of matrix coefficients in any 
selection is linearly independent, and furthermore, any function / e ©(G) may be 
written uniquely as a finite linear combination of the matrix coefficients. Thus it is 
always possible to construct a linear map lo : 6(G) —>■ C that gives any desired set of 
values to the corresponding coefficient functions; in particular, there exists an lo such 
that Lj{Da) — U{a). We now show that lu is multiplicative and conjugate preserving. 
Applying the linear map lo to both sides of eq. (|2.ip results in the following identity: 

(2.3) u;{D, ^ Ds) ^ CasHDc.^} (B ■ ■ ■ (Bcj{D^J]Cls- 
Substituting the matrices U on the right side and using (|2.2p reveals that 

(2.4) Lo{D„®Ds)^Lu{D^)(E)Lo{Ds). 

Consequently, lo is multiplicative. To prove that lo preserves conjugation, note that 
the equation DaD]^ = I implies that uj{Da)Lo{Dl^) = I for all a. Since the matrices 
U{a) ~ io{Da) arc unitary, we also have Lo{Da)u!{Day = I. Matrix inverses are 
unique, and thus w(£'i) = w(I?J,), showing that u; preserves conjugation. By the 
Tannaka-Krein theorem, there exists a unique g £ G such that U{a) = Lo{Da) = 
Da{g) for ah a. Q 

3. Bispectrum. We use this result to establish sufficient conditions for a func- 
tion to be described uniquely by its bispectrum. It is convenient to establish the 
Fourier transform domain for compact groups. Let {Da}aeg be any selection of irre- 
ducible representations. The Fourier transform of any / in Li{G) is the matrix- valued 
function F, such that for each a G Q, we have 

(3.1) Fia) = / fig)D^ig)Ug. 

Jg 

Here the integral uses the Haar measure dg on G; because G is compact, dg is both 
left and right invariant. 
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We use two important properties of the Fourier transform in what foUows ([51 pp 
73-78]): (i) the Fourier transform of any / G Li{G) determines / uniquely up to a set 
of Haar measure zero; (ii) s{g) = r{xg) for aU g if, and only if, S{a) = R{a)Da{x) for 
all a. 

For / e Li(G), the triple correlation a^j is defined as follows: 



(3.2) a3,/(5i,52) = / fiaT f{99i)f{992)dg. 

JG 

(Compare with eq. (|l.ip V Note that the triple-correlation is invariant under left- 
translation, i.e., if there exists x such that r{g) = s{xg) for all g, then a^ ,. = a^^g. This 
follows directly from the left-invariance of the Haar measure dg. Similarly, we may de- 
fine a right-translation invariant version of eq. (|3.2p by integrating f (g)* f {gig) f {g29) , 
but we will not pursue this minor variation in what follows. 

Because / G Li{G), we have that a^j is a function in Li{G xG). It is known that 
any irreducible representation of G x G is equivalent to a tensor product Dcr (8 Dg, 
where D„ Ds are irreducible representations of G (jJTl pg 45]). Thus the Fourier 
transform of 03,/ with respect to the selection {-DajaeS tti^ function on Q x Q that 
is defined as follows: 

(3.3) A3jia,5)^ [ [ a3j{gi,g2)[D^{gi)^ ®Ds{g2)^dgidg2. 

JG JG 

There exists a convenient formula for computing A^ j. 

Lemma 3.1. For any pair a, 5, let C^rs be the matrix and let ai, . . be the 
indices appearing in eq. i2.1]) . Then 



(3.4) A3j{a, 6) - [F{<j) «. F{S)] C^s [F{ai)^ ® • ■ • ® F{ak)^] 



Proof. Since 03 j is integrable, we use the Fubini theorem to interchange the order 
of integration in the following derivation: 

^3,/(o-,^) = / / a3j(5i,52) [D„{gi)'^ ® D„{g2)'^] dgidg2, 

fig)* f {991) f {992) [D„{gi)^ ®Ds{g2)^] dgdgidg2, 

G JG JG 



G JG 



f{9r / / f{99i)f{992) [DA9iy ® Ds{92y] dg^ dg2 dg. 

G JG JG 

By making a change of variables, we find that the double integral inside simplifies as 
follows: 

/ / f (991) f {992) [DA9i)^ ® Ds{g2)^ dg, dg2 = [F{a) ® F(5)] [D,{g) ® Ds{g)] . 

JG JG 

Upon substituting into the expression for A^j^ we find that 
(3.5) A^j = [F{a) ® F{5)] f /(g)* [D^{g) Ds{g)] dg. 



G 



Upon substituting the tensor product decomposition (|2.ip into the above, we obtain 



(3.6) A3j{<J,S) = [F{a)(E,F{6)]C„s 



f{gr{D^,{g)®---®D^,{g))d9 



G 



6 



R. KAKARALA 



After evaluating the integral, the result (|3.4p follows. □ 

The lemma helps to quickly establish the basic completeness result for the bis- 
pectrum on compact groups. 

Theorem 3.2. Let G he any compact group, and let r in Li{G) be such that its 
Fourier coejficients R(a) are nonsingular for all a £ Q. Then a^^g = a^^r for some 
s G Li{G) if and only if there exists x G such that s{g) — r{xg) for all g. 

Proof. If s{g) = r{xg), then the translation-invariance of the triple correlation 
implies that 03^^ = a-3,s- We now prove the converse. Let s be such that 03^^ — a^y^ 
then A^ r = ^3^^, and by Lemma l3.H we obtain that for a, S that 

[R{a) (E> R{6)] C„5 [R{ai)^ © • • • © i?K)^] = 

(3.7) [S{ct) ® S{5)] C,s [5(ai)^ © • • • © 5K)^] 

Set a = 5 = 1, where 1 is the trivial representation g 1 of G. Both i?(l) and S{1) 
are complex numbers, and the equality above becomes 

(3.8) R{l)R{\)R{\f = S'(1)S'(1)S'(1)*. 

Thus i?(l) = <S'(1). Now set 6 — 1\ for any cr, we have D^r ® Di — Da, and thus 
eq. (|3.7p becomes 

(3.9) [i?(cr) i?(l)] R{ay = [S{cr) S{1)] S{a)K 

By assumption i?(l) = S{1) is a non-zero scalar, and we cancel it from both sides 
to obtain that i?(cr)i?((T)t = 5(ct)5((t)^ for all cr. Such an equality between matrices 
holds if and only if there exists a unitary matrix U{a) such that S{a) — R{a)U{a). 
Substituting for S in (|3.7p yields, upon rearranging terms, 

[R{a) © R(5)] G,s [R{ai)^ © ■ • ■ © R{ak)^] = 
[R{(7) R{5)] [U{(t) U{5)] a* [U{ai)^ © • ■ • © U{ak)^] [R{ai)^ © • • ■ © R{ak)^] 

We cancel the nonsingular matrices R{(t) from both sides and rearrange the remaining 
terms to obtain the identity 

(3.10) U{a) © U{5) = Gas [U{ai) © • • • © U{ak)] 

Since the identity above holds for all cr, (5, Theorem 12.21 implies that there exists 
X e G such that U{a) = D„{x) for all cr. Thus S{a) = R{(j)D„{x) for all cr, and the 
translation property of the Fourier transform now implies that s{g) = r{xg) for all g. 
□ 

The hypothesis that all coefficients R{cr) are nonsingular is satisfied generically, in 
the sense that almost every n x n matrix is nonsingular with respect to the Lebesgue 
measure on the set of n x n matrices. Nevertheless, it is desireable to weaken the 
hypothesis, to include for example functions on G that are invariant under the trans- 
lations of a normal subgroup N of G. We prove a result for this case. 

We review some facts concerning group representations and normal subgroups 
(0 Pg 64]). Let iV be a closed normal subgroup of G. Any irreducible representation 
D of the quotient group G/N extends to an irreducible representation of G by 
composition: D — D o tt, where tt is the canonical coset map tt : G ^ G/N. The 
converse is also true: any representation D of G such that D{n) = I for all n G is 



BISPECTRUM ON COMPACT GROUPS 



7 



of the form D = D o tt ior some representation D of G/N. Moreover, letting [G/N) 
represent the dual object of the group G/N, the set 

(3.11) g[N] = {D = boTT,De(GjN^}, 

is closed under both conjugation and tensor- product decomposition, i.e., the tensor 
product of any two representations from the set decomposes into irreducible repre- 
sentations that are also contained in the set. Conversely, to each subset A of Q that 
contains 1 and that is closed under both conjugation and tensor-product decompo- 
sition, there corresponds a unique closed and normal subgroup N of G such that 

A = g[N]. 

Now let / be a function in L2{G) that is invariant under N , i.e., f{ng) = f{gn) = 
f{g) for all n G N. li a ^ then the Fourier coefficient matrix F{a) is a zero 

matrix. To prove this, note that the Peter- Weyl theorem implies that the matrix coef- 
ficients dP^{-) from any selection {^ctlo.gjg/^v-, form an orthogonal basis for L2{G/N); 

thus the corresponding functions = d^'' o tt on G form an orthogonal basis for the 
closed subspace in L2(G) of functions invariant under N. Consequenctly, any TV- 
invariant function in L2{G) has zero inner product with the coefficients of Da when 
a ^ GW]- We use those facts to produce a stronger version of Theorem 13.21 

Theorem 3.3. Let r £ L2(G) be such that its Fourier coefficients R satisfy the 
following conditions: 

1. Each R{a) is either zero or nonsingular; 

2. The set of a such that R{a) is non-singular includes 1, and is closed under 
conjugation and tensor product decomposition. 

Then there exists a normal subgroup N of G such that r is N -invariant, and further- 
more r is uniquely determined up to left translation by its bispectrum A^j. 

Proof. As discussed above, the set of a such that R{a) is nonsingular corresponds 
to Q[N] for some normal subgroup N of G. Furthermore, r — ron for a unique function 
r on G/N. We obtain A3 p from A3 ^ by restricting the latter to the arguments {a, S) 
for which R{a) and R{S) are nonsingular. Theorem 13.21 now shows that r is uniquely 
determined up to a left translation by A3 jr, and thus r = r o tt is uniquely determined 
up to a left translation by A3 □ 

Remark. The hypotheses of Thms 13.21 and 13.31 have an interesting interpretation 
in the context of the Tauberian theorems for compact groups. The latter theorems 
determine what functions lie in the span of translates of a single function / in Li[G). 
Edwards ([21 pp 121-125]) describes one such result: If /i, /2 in £i(G) have Fourier 
transforms Fi, F2, such that F2{a) — Fi{a)M{a) for each a, where M{a) is an 
arbitrary matrix whose dimensions match that of Fi{a), then f2 lies in the span of 
left translates of /i, i.e., /2 may be approximated arbitrarily closely in Li by linear 
combinations of left translates of /i . Suppose now that /i satisfies the hypothesis of 
Theorem [321 i-e., Fi{a) is nonsingular for all a. Then the aforementioned Tauberian 
theorem implies that any function /2 G Li(G) lies in the span of translates of /i, i.e., 
the translates of fi span Li (G) . Similarly, if /i satisfies the hypothesis of Theorem l3.3i 
then the translates of /i span the closed subspace of LiiG) that consists of functions 
invariant under some fixed normal subgroup N . As our theorems show, the bispectrum 
of /i identifies exactly which functions are its translates. 

4. Homogeneous spaces. The definition of a homogeneous space is as follows. 
Let G be any topological group and X any topological space. We say that G acts 
(on the right) on X if for each g E G there exists a homeomorphism Tg : X X, 
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such that Ti,{x) — x for the identity e in G, and furthermore, for gi, 52 in G, we have 
''9132 (^) — ''92 (''si(^))- The group G acts transitively on X if for each xi, X2 in X, 
there exists g d G such that Tg{xi) = a;2. The space X is a homogeneous space for G if 
G acts on X transitively and continuously. An important example of a homogeneous 
space is the quotient space of right cosets G\H ~ {Hg : 5 G G} of a closed subgroup 
H in G. In fact, it is a theorem that any locally compact homogeneous space X of a 
separable and locally compact group G can be represented as a quotient space G\H 
for some closed subgroup H of G ([H pg 124]). 

Our goal in this section is to investigate the bispectrum's completeness for func- 
tions on arbitrary homogeneous spaces of compact groups. By the result cited above, 
we lose no generality by focusing on spaces of the form G\H, where G is some com- 
pact group and H some closed subgroup of G. To any function / on G\H there 
corresponds a unique function f on G such that / = / o tt, where tt : G — > G\H is 
the canonical coset map; conversely, to any function f on G that is invariant under 
left _ff-translations, i.e., f{hg) = f{g) for all g S G and h € H, there corresponds a 
unique function / on G\H such that f = f on. Thus we lose no generality by further 
restricting our study of functions on homogeneous spaces to functions on G that are 
left iJ-invariant for some closed subgroup H . 

Our main tool for proving completeness results is the Iwahori-Sugiura duality 
theorem for homogeneous spaces of compact groups [9]. Let G be any compact group, 
{Da}a£g be any selection of irreducible representations, and 0(G) the representative 
algebra of G. For any closed subgroup H of G, let Qh{G) denote the subalgebra 
of 9(G) consisting of functions that are invariant under left 7?-translations. For 
each / e Qh{G), let f{Hg) denote the common value given to elements of the coset 
Hg by /. The algebraic structure of Qh{G) is revealed to a large extent by the 
multiplicative linear functionals uj : Qh{G) C, i.e., algebra homomorphisms of 
Qh{G). The Iwahori-Sugiura theorem characterizes those algebra homomorphisms 
that preserve conjugation. 

Theorem 4.1 (Iwahori-Sugiura). To each algebra homomorphism lu : Qh{G) —>■ 
C that preserves conjugation, there corresponds a unique coset Hg in the quotient 
space G\H such that for all f e Qh{G), 



We describe an equivalent formulation of the Iwahori-Sugiura theorem that is 
analogous to Theorem 12.21 Several preliminary results are required for the new for- 
mulation, with some of the longer proofs being relegated to the appendices. 

Lemma 4.2. Any function f G 0/j(G) can be expressed as a unique finite linear 
combination of the left H -invariant matrix coefficients of a given selection. 

The proof is given in Appendix [Xj 

Let G, H, and {Da}aeg be as before. Let us define a corresponding sequence of 
matrices {Pa}aeg as follows: 



where dh denotes the normalized Haar measure on H. It is easy to show that each 
Pa is a projection, i.e., a self-adjoint matrix such that PqPq, = -Pa (0 Pg 190]). 
Moreover, the projection matrices as defined above inherit some of the tensor product 
properties of the corresponding representations (fSl pg 190]). 



(4.1) 



^(/) = .f{Hg). 



(4.2) 
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Lemma 4.3. Let {Pa}aeG ohove. For each a, S, let C„s he the Clehsch- 

Gordan matrix and ai, . . ., be the indices in the tensor product decomposition in 
eq. ![MIW- Then 

Pa®Ps= C„S [Pai © • • • © PaJ [Pa ® Ps] , 

= [P, ® Ps] C,S [Pai © • ■ • © Pa,] Cls- 

It proves convenient to apply the following similarity transformations to the P matri- 
ces. For each a, let rank(a) denote the rank of Pa, and let /(rank(a)) be the diagonal 
matrix whose first rank(a) diagonal entries (from the upper left) are 1, and the rest 
are 0. Then there exists a unitary matrix U{a) such that ([15, pg 195]): 

(4.3) U{a)PaU{a)^ = /(rank(a)), 

If we apply the same similarity transformation to the representation D^, then it is 
easily seen that 

(4.4) U{a)Da{h)U{aY = [©;r^"^l(/i)] © D^ih), heH. 

In the decomposition above, 1 is the trivial representation of H , and the last term 
is some unitary representation of H that does not contain 1. 
Rather than starting with an arbitrary selection of {Da}ai^g, suppose now that 

we choose one in which each matrix Da{h) is exactly equal to a direct sum where the 

first rank(a) representations that appear in the sum are 1, i.e., 

(4.5) Da{h) = [©;tt'"^l(/i)] © (h), heH. 

We always obtain such a convenient selection (that is what we shall call it henceforth) 
from a given one by applying similarity transformations as described above. For a 
convenient selection, the projection matrices in eq. (j4.2p are simply Pa — /(rank(a)) 
for all a. 

Lemma 4.4. Let {Da}aGQ be a convenient selection and {Pa}aGQ be its projec- 
tions. The nonzero coefficients in the matrices {PaDa} oltc precisely those coefficients 
of the selection that are left H -invariant. 

Proof. Each matrix PaL>a = I{'rank{a))Da has its first rank(a) rows equal to 
those of Da, while the remaining rows are identically zero. Moreover, PaDa{hg) = 
PaDa{g) for all h and g (simply substitute ^^Dadh for Pa and use the translation 
invariance of the Haar measure dh) , and thus the nonzero coefficient functions in each 
PaDa are left i/-invariant. We now show the converse: any left 7?-invariant coefficient 
d^f of Da is one of the nonzero coefficients in PaDa. Left i7- invariance requires that 

dim(Q) 

(4.6) dl^{g) = dPa'{h9)^ J2 <Wdi^i9)- 

i=i 

The linear independence of the coefficients implies that dPf{h) = 1 for all h ii £ — p, 
and d^{h) = for all h ii I ^ p. But the assumption on Da requires that d^{h) — 1 
on H only if p < rank(Q!), and thus any left i7- invariant coefficient d^^ must appear 
in one of the first rank(Q!) rows of PaDa. □ 

Since the left invariant coefficients are a basis for 0^f(G), any linear map 
Lu : QniG) ^ C is uniquely determined by the values that it gives to those coefficients. 
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For each matrix P^Da, the map uj produces a corresponding matrix u}{PaDa). We 
now determine conditions in terms of the matrices lo{PD) under which w is not only 
Unear but also multiplicative and conjugate-preserving. In the following, we use the 
standard inner product < CijC2 >= C1C2 for complex- valued row vectors Ci, C2, and 
the standard norm ||C|| = (< CiC >)'^- 

Theorem 4.5. Let {Da}aeQ be a convenient selection and {Pa}aeg be its projec- 
tions. Any linear map lu : Qh{G) — > C is both multiplicative and conjugate-preserving 
if and only if the following two conditions hold for all a, S, a in Q : 

(4.7) Lu{PaDa)®U0{PaDa) = [Pa ^ PsjCaS [uj{Pc.,Da.,) ® ■ ■ ■ ® uj{Pa.,Dc,,)]Cls; 

(4.8) u;{PMcj{Pc.Dc,)^ = 

In eq. 7[ ), the matrix CaS and the indices ai,. . . ,ak are as in eg. h2.1\) . 
The proof is given in Appendix [B) 

Let / be a function in ii(G) such that f{hg) = f{g) for all ft, in a given closed 
subgroup H of G. The translation property of the Fourier transform ensures that 
each Fourier coefficient F{a) satisfies the identity F{a) — F{a)Da{h) for all h in H. 
Integrating over h, we find that F{a) = F{a)Pa for all a. We say that each Fourier 
coefficient F{a) is of maximal H-rank if the rank of F{a) equals the rank of Pa- We 
now show that if / is any left _ff-invariant function whose Fourier coefficients F all 
have maximal rank, then / is uniquely determined by its bispectrum ^3 j up to a 
left translation. The proof of our assertion uses the standard notation from linear 
algebra [IS]. For each matrix A, let image(A) and ker(A) denote respectively the 
image and kernel of A. For each a £ Q, let Ha denote the Hilbert space on which the 
corresponding representations act. 

Theorem 4.6. Let G be any compact group, and let H be any closed subgroup of 
G. Let r S Li{G) be invariant under left H -translations. If the Fourier coefficients 
{R{a)}aeQ 0,1^ have maximal H-rank, then 03, r = 03,3 for some s £ Li(G) if and only 
if there exists x G G such that s{g) = r{xg) for all g. 

The proof is given in Appendix [C] 

In the theorem above, we did not require that the function s also be left H- 
invariant. (Equality of bispectra may hold regardless of whether both functions are 
_ff-invariant.) Suppose now that two left _ff-invariant functions r, s are such that 
both have maximal _ff-rank coefficients and both have exactly the same bispectrum. 
The theorem just proved demonstrates that under those conditions, there exits x £ G 
such that s{g) — r{xg) for all g. Yet the element x cannot be arbitrary, for s is left 
i?-invariant, and thus s{hg) = s(g), implying that r{xhg) — r[xg) for all ft, S and 
g £ G. But since r is also left 7J-invariant, we must have r{xg) — r{hxg), and thus 
r{xhg) — r{hxg) for all g and ft. The last identity is always satisfied if x lies in the 
normalizer of H in G, which is the subgroup Nh of G defined as follows: 

(4.9) Nh ^ {x eG -.xH = Hx}. 

(The normalizer of H is the largest subgroup Nh of G such that G itself is a normal 
subgroup of Nh-) In fact, we show that x must lie in Nh in the following theorem. 

Theorem 4.7. Letr, s inLi{G) be two left H -invariant functions whose Fourier 
coefficients R(a) and S{a) both have maximal H-rank for all a. Then a^ r = 0,3, s if 
and only if s{g) — r{xg) for some x G Nh- 

Proof. The "if" assertion is shown above, so we prove the "only if" part. Suppose 
that 03^^ = 0,3,3, and that r, s, both have maximal ff-rank coefficients. Under those 
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conditions, Theorem 14.61 shows that there exits x ^ G such that r{g) = s(xg) for 
all g. Then R{a) = S{a)Da{x) for all a £ Q. Furthermore, the left invariance of 
r implies that R{a) — R{o)Pa for each a, Thus S{a)Da{x) = S{a)Da{x)Pa for 
each a, and combining that with the identity S{a) = S{a)Pa yields S{a)PaDa{x) = 
S{a)PaDa{x)Pa, and thus S{a) [PaDa{x) — PaDa{x)Pa] — 0. By the maximal H- 
rank hypothesis, we obtain that 

(4.10) PaDaix) = PaDa{x)Pa. 

Since P^ = /(rank(Q!)) for a convenient selection, the element x satisfies the above 
equality if and only if the unitary matrix Da{x) is the direct sum of two smaller 
unitary matrices, the first with dimensions rank(Q:) x rank(a) and the second with 
dimensions {n — rank(a)) x (n — rank(a)). For such an x, it follows for any h £ H 
that 

(4.11) PaDa{x)'^Da{h)Da{x) = PaDaix'^hx) = P„ . 

But we now see by Lemma [4.41 that PaDa{x~^hx) — Pa if and only if x^^hx e H. 
The last inclusion holds for all h £ H, and thus x^^Hx — H, or equivalently, x £ Nh- 
□ 

In the interesting special case when G is the group 5*0(3) and H is the subgroup 
of rotations that fix the z-axis, we have that Nh = H. In that case, if r is any 
left 7?-invariant function with maximal i?-rank coefficients, then there are no other 
left iJ-invariant functions with the same bispectrum besides r itself. However, that 
does not mean that the bispectrum uniquely determines r: any function s such that 
s{g) = r{xg) on G has the same bispectrum, although s is not necessarily if-invariant. 

If G = SO (3) and H as above, then the maximal iJ-rank condition is easy to 
satisfy. Here it is well-known that rank(Pct) = 1 for all a £ G ([IS])- Thus an 
arbitrary left i?-invariant function r has maximal _ff-rank coefficients if for all a, the 
matrix R{a) contains at least one nonzero coefficient. That is evidently true if any 
noise is present in measuring r. 

5. Reconstruction algorithms. The completeness theory for arbitrary com- 
pact groups in the preceeding sections can be refined further for the special case when 
the group is SU{2), which is the group of all 2 x 2 unitary matrices with determinant 
-|-1. The group SU(2) arises frequently in applications because it is a double-covering 
of the rotation group SO (3), and in many problems it is more convenient to model 
three-dimensional rotations by elements of SU (2) rather than the corresponding ele- 
ments of S'0(3)-one reason is that the addition of rotations is much simpler for SU (2) 
(Cayley-Klein parameters) than for S0{3) (Euler parameterization). The represen- 
tation theory of SU (2) is known in extensive detail, and we take advantage of the 
special properties of 5J7(2)'s irreducible representations to analyze the bispectrum of 
bandlimited functions. The latter are functions whose Fourier coefficients are identi- 
cally zero except for a finite set of indices. One reason why bandlimited functions are 
important is that any L2-function can be approximated as closely as desired in the 
i2-nietric by a bandlimited function. 

The irreducible representations of SU{2) have several properties that simplify 
our analysis of the bandlimited case {\lBl Chapter 2]). First, there exists one and 
only one irreducible representation (up to equivalence) in each dimension. It is thus 
possible to index the set of all irreducible representations (modulo equivalence) by 
the nonnegative integers, in such a way that for each £ > 0, the representation Di has 
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dimension i + 1. With that indexing, Dq is the trivial representation 5 ^ 1, and Di 
is the self represntation g t-^ g. Furthermore, for any nonnegative integers p, q, the 
tensor product Dp ® Dq reduces explicitly as follows: 

(5.1) Dp®Dq = Cpq [Dp+q © Dp+q^2 ® © ■ ' ■ ® i^lp-g]] ■ 

The unitary matrix Cpq above is the Clebsch-Gordan matrix for p and q. 

We now review briefly the methods of Fourier analysis on SU{2). Let / be any 
L2 function on SU{2). Its Fourier coefficients are the matrices 

(5.2) F{i) = / f{g)D,{g)^dg- I > 0. 

Jg 

The function / is the limit in L2 of the series 

00 

(5.3) Y.(^+l)Ti[F{e)D,{g)]. 

By Lemma [3. II we know that the bispectrum on SU{2) has the form 

Asjip, q) = [Fip) ® F{q)] Cpq [Fip + g)t © F(p + q - 2)t © • • • © F{\p - q\)^] 0^. 

We now devise an algorithm for recovering any real-valued bandlimited function 

with nonsingular coefficients from its bispectrum. Our algorithm makes use of the 

following facts from matrix theory. First, any positive definite matrix H has a unique 

1 j_ j_ 

"positive square root", i.e., a positive-definite matrix such that H^H^ = H f [T51 

pg 181]). The square root is constructed explicitly by diagonalizing H, i.e., finding a 

unitary matrix U such that H — UDW , where D is the diagonal matrix of eigenvalues 

- 1 1 

(here, all nonnegative), and setting = UD^Jj'', where D2 is the diagonal matrix 

containing the positive square roots of the eigenvalues. Although the positive square 
i_ 

root is unique, there are in fact several possible matrix square roots, each formed 

by setting Hq = UQU\ where U is the unitary matrix reducing H to diagonal form, 
and Q is a diagonal matrix whose entries are either the positive or negative square 
roots of the eigenvalues Ai, . . ., A„ of H: 



(5.4) diag^±VAi,...,±VA„ 

It is easy to see that HqHq = H for any such Q. 

The second fact from matrix theory that we use is that any nonsingular matrix A 

1 

has a unique polar decomposition A = Hj^-U ^ where — {AA'^)^, and J7 is a unitary 
matrix. The polar decomposition is unique in the sense that ii A = H^U = H'_^U' for 
positive definite matrices H^, H'_^, and unitary matrices C/, U' , then = H'^ and 
U = U' . It is easy to see that as chosen above is such that det(i?_|_) = |dctA|. If 
the determinant of A is real, then we may choose a square root H of (^^^) such that 
A = HU, where U is unitary and det(iJ) = det(A). The last observation becomes 
important in our analysis of SO(3) below. 

We require one last fact: the coefficient matrix F{1) of any real- valued function 
/ on SU{2) has nonnegative determinant. To see that, recall from above that Di is 
the self-representation of SU{2), and thus 



(5.5) F{1) - / f{g)D^gUg = / f{g) 



G 



d\\gY -dfigy 

df{g) d\\g) 



dg. 
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By evaluating the matrix coefficients above, and using the assumption of real /, we 
find that det[F(l)] > 0. 

Putting all the facts above together, we obtain the following result. 

Proposition 5.1. Let L > 0, and let / be any real-valued function on SU{2) 
whose Fourier coefficients are such that F(£) is a nonsingular matrix for each £ < L, 
and furthermore, F{i) = if i > L. Then f can be uniquely recovered up to a left 
translation from its bispectrum A^j. 

Proof. Since / is real-valued, it follows that F{Q) is a real number. Equation ^ 
shows that ^3 /(0,0) = ^"(0)^, and thus we obtain F{0) by taking cube roots. By 
assumption, F(0) is nonzero, and thus we obtain from (O that 



(5.6) 



^3j(l,0) 

m 



= F(l)F(l)t 



The matrix on the right hand side above is positive definite. Let F(l) — (^--^p(o)^) ^ 

be the positive square-root as constructed above. In polar form F{1) = F{1)U, and 
thus F(l) = F(l)C/'f. The determinant of F{1) is positive, as is the determinant 
of F(l), and thus det [W] = +1. Consequently, W £ SU{2), and we may write 
F{1) = F{l)Di{x) for X = W . If L = 1, then we are done. Otherwise, the following 
algorithm produces matrices F{2), . . . , F{L), such that F{£) = F{t)Di{x) for the 
same x and for all 2 < ^ < L. Since we know F{1) and yl3j(l, 1), we obtain F(2) 
from the upper-left 3x3 submatrix of the following 4x4 matrix: 



(5.7) 



F(l)-i®F(l)-i ^3,/(l,l)C: 



The reason we use the matrix above is as follows. All terms above are known, and if 
we substitute for F[\) and A3j(l, 1), then we find that 



F{1)-^ ®F{1)-^ 



Cji [D^{x)^ Di{x)^] [F{iy^ F(l)-i] [F(l) ® F(l)] Cn [F(2)t © F(0)] C\^C_ 
Cancelling terms, and using the reduction formula (|5.ip . shows that 



(5.8) 



Cii \F{ir^ ^Fil)-"-] ^3./(l,l)Cii = [D2{x)^ F{2)^] ©F(0). 



The upper left 3x3 submatrix of the right hand side is exactly the matrix 
\F{2)D2{x)\\ and we set its adjoint equal to F{2). Having obtained F{2) in that 
way, we obtain F{£) for any £ > 2 from the upper {£+ \) x [£ + 1) submatrix of the 
following matrix 



(5.9) 



^{1-1)1 



F{£ - l)-i F{l)-^\ A^j{£ - 1, l)C(,_i)i 
The same argument as above shows that 

'f{£ - l)-i ® F{l)-A^ A^j{£ ~ 1, l)C(,_i)i = [D,{x)^F{£)^] © F{£ - 2)t. 



^(1-1)1 



We set F{£) equal to the adjoint of the upper (^ + 1) x {£+ 1) submatrix of the right 
hand side, and thus obtain that F{£) — F{£)Di{x). On doing so for all £ < L, the 
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function / on SU{2) obtained by Fourier series expansion with the coefficients F{0), 
F(l), . . ., F{L) is such that f{g) = f{xg) for all g. □ 

It is easy to prove the same completeness result for bandlimited functions on 
SO{3). We describe the few differences that exist, drawing on standard facts about 
representations of SO{3) (18, Chap II]). First, the irreducible representations of 
S0{3) occur only in odd dimensions, and there is exactly one representation (modulo 
equivalence) in each odd dimension. Thus we may list any selection of irreducible 
representations as {-D^j^o' where for each £, the representation D( has dimension 
{2£ + 1). In that indexing, Dq is the trivial representation, and Di is equivalent to 
the self-representation g i-^ g oi SO{3), i.e., Di{g) = UgU^ for some unitary matrix 
U. (Recall that 5*0(3) is the set of all real-valued 3x3 orthogonal matrices with 
determinant +1.) For each n, to, the tensor-product D„ ® reduces explicitly as 
follows: 

(5.10) Dn ® D,n = Cn„i [Dn+m ® -Dn+m-1 © ■ ' ■ © ^nm' 

With the formula above, it is easy to see that the recursive algorithm given in the 
proof of Proposition ! 5 . 1 1 generalizes to recover all real- valued bandlimited functions on 
SO{3). To initialize the algorithm, we require an estimate F{1) of the first coefficient 
F{1) from the data i^(l)F(l)'t', such that F(l) = F{l)Di{g) for some element g of 
S0{3). Assuming that F{1) is nonsingular, we obtain the estimate as follows. The 
representation Di is such that Di{g) = UgU\ where U is fixed as g varies in 5(9(3). 
Thus 



^^(1) = / /(.g)Di(g)tdg, 

JG 



u 



G 

' I{g)g^ 

G 



Let Fg(l) denote the matrix that results by evaluating the integral in brackets. Since 
/ is real- valued, and every matrix g has real coefficients, the matrix J^s(l) has only real 
coefficients. Thus the determinant of F{1) — UFs{l)U^ is a real number. Assume for 
the moment that det [F{1)] = det [Fs(l)] > 0. Let F{1) and Fs{l) denote respectively 
the (unique) positive square roots of -F(l}F(l)t and F,{l)F,{l)'i . Since F(l)i^(l)t = 
UFs{l)F4l)^U\ it is easily seen that F{1) = UFs{l)W. Now consider the polar 
decomposition Fs{l) = HV, where H is positive definite and V is unitary. Recall from 

the earlier discussion for SU{2) that H = (^,(1)^^(1)+)^, and thus H = F,{1). Since 
Fs{l) is real- valued, V must be real- valued orthogonal matrix. Matching determinants 
on both sides of the equation Fs{l) = ^,(1)^ reveals that det[y] — +1, and thus 
V — g, for some g e S'0(3). Substitution reveals that 

(5.11) F(l) = UFs{l)U^ = UFs{l)gU^ = UF,{l)U^UgU^ = F{l)D,{g). 

The assumption that det[F(l)] > is not critical. We use it only to obtain 
that det[V] = +1, where V = F^^l)-'^ Fsil). Instead of selecting F(l) to be the 
positive definite square root of F{l)F{iy , we may choose F{1) to be any square 
root such that det[_F(l)] — det[F(l)], e.g., by multiplying the top row of the positive 
definite square root matrix by —1 if necessary. We do not know det[_F(l)] a priori, 
but if we store it as "side information" along with the bispectrum, then we obtain 
a complete rotation-invariant description for any real-valued bandlimited function on 
S'0(3). Note that det[i^(l)] remains invariant under translation on SO{3), i.e., if 
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f{g) = s{hg), then F{1) = S{l)Di{h), but since det[Di{h)] = +1, we obtain that 
det[i^(l)] = det[S'(l)]. To sum up, any real-valued bandlimited function / on 5*0(3), 
whose coefficient matrices are all nonsingular up to the bandlimit, can be recovered 
completely — up to a single translation on 50(3) — if both its bispectrum and the value 
of det[F(l)] is known, and the algorithm described above is used. 

6. Applications. As mentioned in the introduction, the invariance and com- 
pleteness properties of the bispectrum lend themselves to applications in pattern 
matching problems. One particular application is described here. R. Kondor [14j 
demonstrates how translation- and rotation- invariant matching of hand-written char- 
acters is accomplished with bispectral invariants. To do so, Kondor notes that, for 
practical purposes, the characters themselves may defined as intensity-valued func- 
tions on a compact patch on of radius 1. A transformation may be constructed 
that maps the planar patch to the upper hemisphere of the sphere 5"^ as follows: 

(6.1) uj : (r, (j)Tg_2) i—>- [9, (I)g2) with r = sm{d); (j)^2 — (j) gi . 

The subscripts denote the domain of the angle involved, whether it be the plane 
or the sphere S'^ . Kondor shows that rigid body motions on the patch, each of which 
consists of a rotation by an angle a and a translation by a vector T = {tx,ty) with 
\\T\\ < 1, may be mapped to 3-D rotations through the use of Euler angles {O.^'^ip) 
as follows 

(6.2) a = ip; = sin 9 cos (j)', ty = sin 9 sin (f). 

This mapping produces a local isomorphism between planar rigid motions and spher- 
ical rotations. 

By using the transformation (j6.ip . every intensity function defined on the planar 
patch may be converted to a function on 5^. The problem of finding rigid- motion 
invariants on now becomes one of finding rotation invariants on the sphere 5^. 
Since the sphere is a homogeneous space for 50(3), every function / on 5^ may in 
turn be lifted to a function / on 50(3) using the "north-pole" mapping: if z = [0, 0, 1], 
then f{R) = f{Rz) for every R G 50(3). We may now construct the bispectrum of / 
from eq. p.4p using the Fourier transform on 50(3), which may be calculated using 
spherical harmonic basis functions. Kondor calculates bispectral invariants in this 
way, and shows, in an experiment using 1000 hand-written characters from a standard 
dataset, that the invariants perform well in matching over arbitrary orientations and 
starting positions of the characters. 

A second application of the bispectrum occurs in astrophysical models of primor- 
dial fluctuations, as mentioned in the introduction. Cosmic inflation ^ predicts a 
Gaussian pattern of temperature anisotropics in the cosmic microwave background 
radiation (CBR). The CBR anisotropy is a function defined on 5^, and therefore we 
may calculate its bispectrum using eq. (|3.4p . If the anisotropy is Gaussian, then the 
expected value of the angular bispectrum is zero. However, as X. Luo [IS] shows, 
the stochastic nature of anisotropics means that cosmic variance makes it difficult to 
extract non-Gaussian structure from CBR data. In that paper, as in much of the 
physics literature, expressions of the bispectrum follow the "summation notation" , 
which implicitly focuses attention at the level of individual elements of bispectral ma- 
trices. It is hoped that the approach of this paper, including the matrix form derived 
in (j3.4p . proves useful in allowing insight into higher level properties, such as matrix 
rank, decomposition, and completeness. 
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Both of the apphcations mentioned apply the bispectrum to functions on the 
sphere. Healy et al [IDj describe a fast "divide-and-conquer" discrete Legendre trans- 
form, which leads to an "FFT" on S'^. They show how this transform leads to efficient 
computation of bispectrum on the sphere. 

7. Summary and future directions. This paper derives completeness proper- 
ties of the bispectrum for functions defined on compact groups and their homogeneous 
spaces. A matrix form of the bispectrum is derived, and it is shown that every function 
with nonsingular coefficients is completely determined, up to a group translation, by 
its bispectrum. A reconstruction algorithm for functions defined on the groups SU (2) 
and 5*0(3) is described. 

Results similar to those in this paper may be established for non-compact, non- 
commutative groups [H] . Those results rely on the duality theorem of Tatsuuma. The 
Tannaka-Krein duality theorem, which is central to this paper, has been extended to 
compact groupoids [1]. It would be interesting to see if a corresponding bispectral 
theory may be constructed there. 

Appendix A. Proof of Lemma 14.21 Since 9^ C 6(G), each / e Qh{G) 
is a unique finite linear combination of matrix coefficients d^'' (not necessarily H- 
invariant): 

(A.i) fig) - Cdl'ig). 

a,p,q 

We set f{hg) — f{g) to obtain 

(A.2) Yl ^'dl'(9) = E ^I'dl'^hg). 

a,p,q Oi.p.q 

The multiplication rule D^Qig) — Da{h)Da{g) for representation matrices reveals 
that 

dim(a) 

(A.3) d^a'ihg)^ J2 dl'^di^g), 

1=1 

and thus 

(A.4) cS'dS' (5) = E '^''d'' (''X' (9) + E E < (^)^« (5) • 

a,p,q oi,p,q a,p,q 

The linear independence of the matrix coefficients implies that in the equation above, 
dvP{K) = 1 and dv^(h) = if ^ ^ p. Thus d^^(hg) = d^^{g) for every coefficient 
function in eq. (|A.1[) . 

Appendix B. Proof of Theorem 14.51 If uj preserves multiplication and 
complex-conjugation, then the Iwahori-Sugiura theorem shows that there exists a 
unique coset Hg such that uj^PaDa) — PaDa{Hg) for all a. From this, equations 
(|4.7p and (|4.8p follow immediately. Suppose now that ui is some linear map that also 
satisfies eq. (|4.7p . Applying Lemma to both sides of the tensor product decompo- 
sition in eq. (|2.ip yields 

(B.l) (P^D,) ® (PsDs) - [Pa ® Ps] CaS [{Po^^D^J © • ■ • ® (Pa,i?aJ] Clg. 
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Now apply lo to both sides to obtain 

uj (iP,D„) ® (PsDs)) = [Pa Ps] C,5 [u{P„,D^,)® • ■ • ® Lo{Pc.,D^^)] Cl^. 

Because eq. (|4.7I) holds, we obtain that for all cr, 5, 

(B.2) UJ {{PaD,) {PsDs)) = uj{P,Da) ® lu{PsDs). 

Thus is multiplicative. Suppose now that the linear and multiplicative map lj also 
satisfies eq. (|4.8I) . Applying cj to both sides of the identity {PaD„){PaD„Y = P^ 
yields 

(B.3) uj{P^D^)w {{PM^) = P,. 

We show that uj ((PaZ^a)^) = ^^(Pai^Q)^, proving that cj preserves conjugation. Let 
(a be any nonzero row of the matrix PaDa- We establish the following three equalities: 

(B.4) <c^(C„),c^(Ca)> = l, 

(B.5) <c.(Cc.),c.(C)*> = l, 

(B.6) <ujiC),u;iC)> = l. 

The first equality (jB.4p follows from (|4.8p (recall that we are working with a convenient 
selection, for which Pa = /(rank(a))). The second is derived from eq. (|B.3p . The 
final equality requires more work, but is a straightforward consequence of (j4.8[) and 
the linearity of lo. We give its proof later, but for now assume that it is true. The 
three inequalities above imply that 

(B.7) <^(c.),^(cr >=iic.(c.)iiii^(crii. 

The Cauchy-Schwartz inequality shows that the identity above holds if and only if 
uj{(a) = cw(C)*, and from eq. (|B.5p we see that c = 1. Thus w(Cq) = '^(C)*- Since 
the preceeding argument applies to any nonzero-row of any matrix PaDa, it follows 
that CO preserves conjugation. 

Now to prove (|B.6p . For any representation Da in our selection, the conjugate 
representation D* is also irreducible, and there are two cases: (i) -D* = AaDaA}^ for 
some unitary matrix Aa; (ii) -D* = AjiDpA^^ where (3 ^ a. Assume that the first case 
is true. It is easy to show that any matrix Aa expressing the equivalence of conjugate 
representations is symmetric, and thus = A*a (O pg 15]). Furthermore, only the 
first rank(Q;) rows of Da are i?- invariant, and that must also be true of the matrix 
Da = AaDaA*a. Thus Aa transforms the first rank(a) rows among themselves, which 
means that Aa must have the block-diagonal form Aa = Aa^i ® Aa^2i where Aa^i is a 
symmetric unitary matrix with dimensions rank(a) x rank(a). Thus PaAa = AaPa- 
Putting those facts together, we obtain the following identity by virtue of w's linearity: 

(B.8) ^{Pc.D*a) = Oj{PaAaDaA*a) = AaLo{PaD a) A*a. 

By using the identity above and eq. (|4.8p . we find that u}{PaD*a)u}{PaDa)'^ — Pa- 
Noting that Pa — /(rank(Q;)), the previous equality for matrices implies eq. (|B.6[) for 
their nonzero rows. Case (ii) is similar. 

Appendix C. Proof of Theorem [4T6l 
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If r and s are left translates of each other, then 03 — a^,s, as follows from the 
definition of triple correlation and the left invariance of Haar measure. Now suppose 
that as.r — 03. s- Lemma |3 . 1 1 shows that for all a, S, 

(C.l) [R{a) ® R{S)] C,5 [R{.ai)^ © • • • © Cls= 

[S{a) ® S{5)] Ca, [5(ai)t © • • • e 5(afc)t] . 

If we set a = b = \ and apply the same argument used in the proof of Theorem 13.21 
we obtain i?(l) — S'(l). The maximal _ff-rank assumption implies that the scalar 
-R(l) = >5'(1) is nonzero. Now set (5 = 1 in (|C?2)) above. Cancelling R{\) = S{\) from 
both sides shows that 

(C.2) R{a)R{a)'^ = S{a)S{a)'^ . 

Hence, we have for each a that S((j) ~ R{a)U{a) for some unitary matrix U{a). 
Substituting into eq. (|C.2|) reveals that 

[R{a) R{6)] C^s [Riai]^ © • ■ • ® Riak^] = 
[Ria) ® R{S)] C„s [[/(a)t © ... © U(ak)^] [i?(ai)t © . . . ® R{ak)^] Clg. 

We cancel C^^ from both sides. The identity R(a) — R(a)Pa implies that R{ay = 
PaR{a)\ and thus image(i?(Q:)t) c Pai'Ha)- But the rank of equals that of 

Pa, and thus image(i?(Q:)t) = Pa{'Ha)- The last identity implies that 

[R{cj) ® R{5)] C,s [Pai © • ■ • © Pc,] = 
[R{a) ® R{5)] [U{a) ® U{5)] C^s © ... © C/K)t] © . . . © p„J . 

Multiplying both sides from the right by Clg[Po- © Ps] and using Lemma [4.31 reveals 
that 

[R{a) © R{d)] [Pa © Ps] = 
[Ria) © RiS)] [Uia) © {/(J)] C^s © • ■ • © Uiak^] [Pa © /^i] . 

We substitute R{a) — R{a)Pa and R{6) = R{d)Pcr into the leftmost tensor product 
term on the right hand side of the equation above and simplify, to obtain 

[R{a) © R{d)] [Pa © Ps] = 
[Ria) © R{6)] [Pa © Ps] [Uia) © C/((5)] C^s [t/(ai)^ © • ■ • © C/(afc)t] [P, © F5] . 

For each a, the identity Ria)Pa ~ Ria) together with the assumption that i?(a) 
has maximal iJ-rank imply that i?(a) is one-to-one on Pai'Ha)- Thus the equation 
above implies that 

Pa®Ps = [Pa ® Ps] [Uia) © Ui5)] Cas [C/(ai)t «... © Uiak)^] [P^ © Ps] 

The matrix in between the two orthogonal projections on the right hand side is uni- 
tary; it is easily seen that for any unitary matrix U and any orthogonal projection P, 
the matrix equation P — PU P holds only if UP = P. Thus we obtain 

(C.3) Pa®Ps^ [Uia) © Ui5)] Cas [Uiai]^ © ... © Uiau)^] [P^ © Ps] 
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Rearranging terms, we obtain the following identity: 

[f/(a)t C/((5)t] [P„ ® Ps] = C,s [U{ai)^ © • • • © U(ak)^] C^JP, ® Ps]. 

Substituting from Lemma 14.31 in the right hand side, and subsequently taking the 
matrix adjoint of both sides, reveals that 

[(P,C/(a)) {PsU{5))] = [P, ® Ps] Ca& [(P„,C/(ai)) • • • © {P^JJ{ak))\ Cl^. 

Theorem 14.51 together with the Iwahori-Sugiura Theorem, shows that the identity 
above holds if and only if there exists a coset Hx such that PaU{a) — PaDa(Hx) for 
all a G tj. Thus for each a we have the string of identities 

(C.4) S{a) ^ R{a)U{a) = R{a)PaU{a) = R{a)Pa,Da{Hx) = R{a)Da{x). 

The translation property of the Fourier transform now shows that s{g) = r{xg) for 
aU g. 
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